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The gauge fixing procedure for N = 1 supersymmetric Yang-Mills theory (SYM) is 
proposed in the context of the stochastic quantization method (SQM). The stochastic gauge 
fixing, which was formulated by Zwanziger for Yang-Mills theory, is extended to SYM4 in the 
superfield formalism by introducing a chiral and an anti-chiral superfield as the gauge fixing 
functions. It is shown that SQM with the stochastic gauge fixing reproduces the probability 
distribution of SYM4, defined by the Faddeev-Popov prescription, in the equilibrium limit 
with an appropriate choice of the stochastic gauge fixing functions. We also show that the 
• BRST symmetry of the corresponding stochastic action and the power counting argument 

in the superfield formalism ensure the renormalizability of SYM4 in this context. 

§1. Introduction 

00 : 

' In a previous paper/) we applied the stochastic quantization method (SQM)^) 

■ to the M = 1 supersymmetric Yang-Mills theory (SYM). In 4 dimensions, the appli- 

ed ! cation of SQM to SYM4 is formulated in terms of the superfield formalism. Based 

' on the Ito stochastic calculus, the underlying stochastic process, which is described 

by the superfield Langevin equation and also by the corresponding Fokker-Planck 
, equation, preserves the global supersymmetry as well as the local gauge symmetry. 

^j' We have postponed the detailed explanation of how SQM recovers the conventional 

^ ■ (path-integral) probability distribution constructed by the Faddeev-Popov prescrip- 

tion. The main purpose of this paper is to present a proof of the equivalence of the 
SQM approach and the BRST invariant path-integral approach, and to demonstrate 
the perturbative self-consistency of the gauge fixing procedure in SQM for SYM4. 
5^ I In the SQM approach, as far as the expectation values of gauge invariant observ- 

ables (such as Wilson loops) are concerned, it is not necessary to introduce an explicit 
gauge fixing procedure. By contrast, for gauge variant quantities, we need a gauge 
fixing procedure to introduce a drift force for the convergence of the longitudinal 
degrees of freedom. Such a procedure, the so-called stochastic gauge fixing proce- 
dure, was first formulated for Yang- Mills theory (YM),^) and its equivalence to the 
Faddeev-Popov prescription in the path-integral approach was shown. It was also 
clarified how the extra gauge fixing term simulates the Faddeev-Popov ghost loops in 
the perturbative sense. ^^'^^ The renormalization procedures in various models have 
been investigated in two different schemes. One is based on the so-called stochastic 
action. ®)"^^) The other respects the gauge covariantly regularized Schwinger-Dyson 
equations. For non-super symmetric gauge theories, in particular, renormalizabil- 
ity in SQM was established in terms of the Ward-Takahashi identities, which are 
derived from the BRST symmetry of the stochastic action. The BRST invariance 
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in this context has also been studied in relation to the five-dimensional* ^ local gauge 
invariance of YM4,^'^)'-'^^) first-class constraint systems including gravity/^)' and 
a topological field theory approach to SQM.^'^^ 

In this note, we formulate the stochastic gauge fixing procedure for SYM4 in the 
superfield formalism. It is defined by introducing a chiral and an anti-chiral super- 
field as the gauge fixing functions. It is shown that SQM applied to SYM4 with the 
stochastic gauge fixing procedure reproduces the standard Faddccv-Popov probabil- 
ity distribution under the appropriate choice of the gauge fixing functions, while the 
expectation values of local gauge invariant observables are independent of the gauge 
fixing procedure, in particular, of the specific choice of the gauge fixing functions. 
The stochastic gauge fixing procedure enables us to carry out the perturbative analy- 
sis in terms of the stochastic action principle in the superfield formalism. The power 
counting argument indicates the renormalizability of this approach. For SYM4 in 
the superfield formalism, the renormalizability was investigated in the conventional 
BRST invariant path-integral approach. -"^^^ We show the renormalizability of the 
BRST invariant stochastic action principle by applying the standard procedure in 
terms of the Ward-Takahashi identity for the 1-P-I vertices. 

This paper is organized as follows. In §2, we recapitulate the SQM approach for 
SYM4 in terms of the superfield Langevin equation. In §3, we introduce the stochastic 
gauge fixing procedure. Section 4 presents a formal proof of the equivalence of the 
stochastic gauge fixing procedure and the conventional Faddeev-Popov prescription. 
In §5, we introduce the stochastic action in the BRST invariant framework. We show 
that there exists an extended BRST symmetry in a fivc-dimcnsional sense and that 
it is truncated to a nilpotent reduced BRST symmetry by integrating out auxiliary 
fields. In §6, we give the power counting argument and sketch the renormalizability 
of the BRST invariant stochastic action principle. Appendix A is devoted to explain 
the conventions to make the paper self-contained. In Appendix B, we give detailed 
discussion to complete the proof of the renormalizability of the BRST invariant 
stochastic action principle in terms of the Ward-Takahashi identities. 

§2. Superfield Langevin equation for SYM4 

In the SQM approach for SYM4, the superfield Langevin equation and the cor- 
responding Fokker-Planck equation are formulated in terms of the Ito stochastic 
calculus by exploiting the following similarity of SYM4 in the superfield formalism 
to lattice gauge theories:^) 

dynamical field : Ufj, G U{N), link variable 

^ U = e^^^ , vector superfield V, V € su{N) . 
local gauge transformation : ei^{^)u^e~i-^i^+p) 

^ e-'^'^^Ue's^, chiral supcrfields DA = DA^ = . 
diflFerential operator : Ea, (left) Lie derivative, EaU^ = taU^ 



There, the stochastic time t plays the role of the 5-th dimensional coordinate. 



'short' title 



3 



^ £ai analogue of the (left) Lie derivative . 
integration measure : dU^, Haar measure on U{N) 
<-> VGW, analogue of the Haar measure . 

In lattice gauge theories, the stochastic process that describes the time evolution 
of link variables is well-established.^'')'^^) Although U{z) = e^^^^^^ is not a group 
element, the correspondence can be realized by introducing a differential operator 
^a^"* {z) ( '^1^'' {z) ) and a local gauge invariant measure ^/G'DV , which are analogues 
of the left (right) Lie derivative and the Haar measure defined on a group manifold, 
respectively. This correspondence is useful for the proof of the equivalence of the 
stochastic gauge fixing procedure and the conventional Faddeev-Popov prescription 
given in §4. It also simplifies the application of the standard procedure for the 
renormalizability in terms of the Ward-Takahashi identities to the BRST invariant 
stochastic ation principle. 

The superfield Langevin equation for SU{N) SYM4 is obtained by regarding 
e^^^ as a fundamental variable: 

— (Z\e29^)e-2^'^(r, z) = -fiAT—le^^^^VW + VW) + Aw(t, z) , 

{Awij{T, z)AWmn{r, z'))Awr = '2kpAT(dinSjm - j^Sijdmn^ 6^{z - z') . (2-1) 

Here, {.■.)Awr represents that the expectation value is evaluated by means of the noise 
correlation at the stochastic time r. /3 is a scaling factor for the stochastic time which 
does not affect to the equilibrium distribution. The conventions employed here are 
explained in Appendix A. The Fokker-Planck equation is deduced from the Langevin 
equation (2-1): 

|:P(r,e2^^) = /3^ yd8^TV^(^)(2)(^(^)(z) + (^(^)(z)5))p(T,e2s^) . (2-2) 

Here, P{t, e^^^) is a scalar probability density. Therefore, in a formal sense, 
e"^ is a stationary solution of the Fokker-Planck equation. However, to assert the 
uniqueness of the equilibrium limit, we must introduce a regularization procedure 
to reduce the number of degrees of freedom to a finite value. It is also necessary 
to introduce the stochastic gauge fixing procedure for the convergence of the gauge 
degrees of freedom. 

The Langevin equation and the noise correlation in (2-1) are covariant under the 
stochastic-time independent local gauge transformation e^^^ g-«s-^''^g23Vg«9^^ fg,. 
which A and are stochastic-time independent, provided that the noise superfield 
is transformed as 

Aw{t, z) e~'^^^W^ Aw{t, z) . (2-3) 

This means that Aw is not a vector superfield. The introduction of the complex 
noise superfield Aw^ ^ Aw is the price we have paid for the local gauge covariance 
of the superfield Langevin equation (2-1) for e^^^. 
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The superfield Langevin equation for the vector superfield V is deduced from 
(2-1) as 

AV{r, z) = -(3Ar{^-^^VW + ^^^^m) + A^S{r, z) . (2-4) 

Here, AyjE is the collective noise superfield defined as 

. 2gLv 
A^^{t,z)= ^^^^^ _^ Aw 

= ^|coth(^Ly) — Iwh^gLy) ^Lv{Aw + Aw^). (2-5) 

To derive the second expression, we have used the hermiticity assignment of the 
noise superfield (A- 10). As is clear from this expression, this collective noise is a 
vector superfield, i.e., A^E^ = A^S. Therefore, the time evolution of the vector 
superfield V preserves its hermiticity = V. We note that the drift force term on 
the r.h.s. of (2-4) is not a simple variation of the action, 6S/6V"', but it includes a 
superfield kernel G°-^5S/5V^ for the local gauge covariance, as explained in (A-9). In 
the weak coupling limit, our formulation is identical to those in Refs.22)— 24) for the 
super symmetric U(l) vector multiplet. 



§3. Stochastic gauge fixing for SYM4 



In the SQM approach, the stochastic gauge fixing was first introduced to solve 
the Gribov ambiguity^^^ in YM4, which exists in a class of covariant gauges in the 
conventional Faddeev-Popov prescription.^)"^) For YM4, the gauge fixed Langevin 
equation is given by 

AAl{r,x) = -^r{^^ - V,.«(yl)} +Zl^^(r,x) , 

{Ar^l{T,x)Ari{T,y)) = 2At5'^^5^X{x - y) , (3-1) 

where Sym = ^ jd^xF'^^F^^ and {V t^vf = d^v'' + /'^^Mj^?;^. v'^iA) is the stochastic 
gauge fixing function. In the Zwanziger gauge fixing, it is chosen to be local: 

v'' = -(d-Ar. (3-2) 
a 

For any functional of the gauge field, the stochastic gauge fixing term appears as the 
generator of the local gauge transformation in the stochastic-time evolution equation 
of their expectation values. Therefore, it does not affect the expectation values of 
local gauge invariant observables. In addition, the significant fact established in 
the perturbative calculations is that the extra self-interaction of the Yang-Mills field 
introduced by the stochastic gauge fixing term in (3-1) simulates the Faddeev-Popov 
ghost contribution in the conventional approach.^)' 

Formally, the equivalence of the stochastic gauge fixing procedure and the Faddeev- 
Popov prescription is shown by specifying the stochastic gauge fixing function in such 
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a way that the Faddeev-Popov distribution, 
Pfp[^] = J VcVcVBe~^° , 

So = Sym + + B'^d ■A'^ + dcf- Vc«} , (3-3) 

satisfies the stationary condition for the probabihty distribution in equihbrium: 

(3.4, 

Here, the stochastic gauge fixing function v^'lx; A] is given by 

+a-'d,{d ■ A)\y) + 25/^'='^(c'=(y)5,c^(y))} , (3-5) 

where the ghost propagator is defined by regarding i^p in (3-3) as the partition 
function. As is clear from the expression of i;'*[.t; A], it is non-local with respect to 
the gauge field, and its explicit form may be determined perturbatively with the 
appropriate renormalization procedure: 

v''[x;A] = -(d-Ay + 0(A'^). (3-6) 
a 

Although the stationary condition (3-4) with (3-2) cannot be solved for an arbitrary 
value of the gauge parameter a, (3-6) indicates that the Zwanziger gauge fixing 
(3-2) corresponds to the non-local gauge fixing in the conventional approach. In 
a perturbative sense, this is the reason why the self-interaction of the gauge field 
in the stochastic gauge fixing term simulates the contribution of the ghost in the 
conventional approach. The gauge fixing function (3-2) ( and also (3-6)) provides 
a restoring force for the longitudinal mode of the gauge field, and therefore, it is 
expected in perturbative calculations that the probability distribution relaxes to the 
Faddeev-Popov distribution in equilibrium. 

In a non-perturbativc sense, however, the Faddeev-Popov distribution is not 
positive definite, because det(9 • V) is not positive due to the presence of Gribov 
copies, and it is not justified that the probability distribution at finite stochastic time, 
P[r; A], uniquely relaxes to the Faddeev-Popov distribution -F^rpfA] in equilibrium. 
The stochastic gauge fixing procedure enables us to introduce a class of covariant 
gauges which are non-holonomic in general. Therefore, if it happens in the stochastic- 
time evolution of the gauge field that the configuration of the gauge field is confined 
to the Gribov region, where the Faddeev-Popov operator —d ■ V{A) is positive, 
the probability distribution P{t) uniquely relaxes to the equilibrium distribution 
Ppp for which the support of the distribution is restricted to the Gribov region. 
Unfortunately, it is unclear whether this is indeed the case or not.^^ In Ref.26), it 
is argued from a slightly different point of view that the solution of (3-4) with (3-2) 
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in the limit of the Landau gauge (a — > +0) can be defined by the Faddeev-Popov 
distribution in the Landau gauge for which the support is restricted to the Gribov 
region f2 = ■ A = 0, —d ■ V{A) > 0}. In the Langevin equation (3-1) with 

(3-2), we may ignore the term 6Sym/SA in the hmit a — > +0. Then we obtain 

A(d ■ A) = Ar-d ■ V(A)(d ■ A) + d ■ At] . (3-7) 

a 

Therefore, if the configuration of the gauge field initially exists inside the Gribov 
region J7, the random walk along the gauge orbit is strongly suppressed in the sin- 
gular limit of the restoring force {a — »• +0), while there exists at least one unstable 
direction along the gauge orbit outside of 1?.^^^ We note that the restoring force 
is non-singular along the transverse direction. In this regard, on the basis of the 
Langevin equation, it is unclear whether the random walk of the gauge field config- 
uration along the transverse direction, i.e. A = d- AA = 0}, is restored 
enough to be confined to f2. In addition, since there are Gribov copies inside the 
Gribov region, in a precise resolution of the Gribov problem, the integration must 
be restricted to the so-called fundamental modular region, which is a proper subset 
of the Gribov region and free of Gribov copies.*-* Although there have been some pro- 
posals for defining a globally valid path-integral for YM4,^^)'^^) the relations between 
these proposals and the stochastic gauge fixing procedure are also open questions. 

In the case of SYM4, it is plausible that we encounter the Gribov problem 
provided that we work with a class of covariant gauges in the Wess-Zumino gauge. 
Although this may cause a problem similar to that which occurs in YM4 , it is beyond 
the scope of this paper to study the Gribov problem in SYM4. In this paper, we do 
not investigate the Gribov problem in terms of our choice of the stochastic gauge 
fixing for SYM4 and restrict our argument to showing the formal equivalence of 
the stochastic gauge fixing procedure and the conventional naive Faddeev-Popov 
prescription by specifying the stochastic gauge fixing function in such a way that 
the Faddeev-Popov distribution of SYM4 satisfies the stationary condition of the 
probability distribution in equilibrium. The main purpose of §3 and §4 is to define 
the stochastic gauge fixing procedure for SYM4 and to derive the supersymmetric 
extension of Eqs.(3-1), (3-2), (3-4) and (3-5). 

In order to define the stochastic gauge fixing procedure for SYM4, we first intro- 
duce the stochastic-time dependent auxiliary superfields A and A'^ , which are chiral 
and anti-chiral, respectively, into the Langevin equation by performing the (inverse) 
local gauge transformation e^^^ g«s-^^g23Vg-«s/l_ Then we redefine the auxiliary 
chiral and anti-chiral superfields as 

0=- ^,<^=^ 7 3-8 

AA 

Here we have denoted A = ——. The redefined auxiliary superfields, 6 and d>, are 
At 

*^ It is also conjectured in Ref.26) that the Landau gauge Faddeev-Popov distribution restricted 

to fi gives the same correlation functions, {A{xi)A{x2)..-A{xn)) , as those defined by the Landau 

gauge Faddeev-Popov distribution restricted to the fundamental modular region. 
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also chiral and anti-chiral, D(p = D(p = 0. Equation (2-1) becomes 

= -pAr^ie^'^^^VW + VW) + Aw . (3-9) 

This Langevin equation is invariant under the stochastic-time dependent local gauge 
transformation defined by 

^igA{T,z) _^ ^igA{T,z) ^igE{T,z) _ (3-10) 

Here, the stochastic-time dependent transformation parameters, E and S\ are chi- 
ral and anti-chiral, respectively. Under this extended local gauge transformation 
( "extended" means stochastic-time dependent)^ the Langevin equation is covariantly 
transformed. In particular, the auxiliary superfields are transformed as 

1 j p-igLs _ 1 . 
-. . —igL^ -. 

^ ^ r ~ + e-'^^^t^ . (3.11) 

P 9 J-'E'f 

The extended local gauge invariance of the Langevin equation (3-9) under (3T0) 
and (3-11) must be broken by the gauge fixing procedure. In order to fix the gauge 
, we specify the gauge fixing functions, (f) and ^, as follows: 

= i-D'^D'^V, ^ = -i-D'^&V . (3-12) 
4 4 

This choice of the gauge fixing functions is determined almost uniquely by dimen- 
sional analysis if we require that the Langevin equation is manifestly supersymmetric, 
invariant under the 4-dimensional Lorentz transformation and that the gauge fixing 
functions are linear with respect to the vector superfield V . In general, we may 
choose 

<\> = i^D'^D'^hiV), 4> = -i^D^D^h{V) , (3-13) 

with an arbitrary function h to specify the local gauge fixing term. The Langevin 
equation (3-9) with the gauge fixing functions (3T2) is a supersymmetric extension of 
the gauge fixed version of the YM Langevin equation. In the component expansion,^) 
(3-9) with (3-12) includes the expression of the YM Langevin equation with the 
Zwanziger gauge fixing term (3-2). Therefore, it provides the stochastic gauge fixing 
that is the supersymmetric extension of the so-called Zwanziger gauge in the YM 
case. We show in the next section that the local gauge fixing functions (3-12) are the 
leading terms of the non-local stochastic gauge fixing functions which are equivalent 
to the manifestly Lorentz covariant local gauge fixing in the conventional Faddeev- 
Popov prescription. The perturbative aspect of this gauge fixing procedure is studied 
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in §5. For example, the superpropagator of the vector superfield, given in (5-15), is 
identical to that derived from the Langevin equation (3-9) with (3-12) in the Feynman 
gauge ^ = 1, excluding the initial condition dependence of the superpropagator in 
solving the Langevin equation. 

Next, we show that the gauge fixing procedure we have introduced does not 
affect the time development of the local gauge invariant observables. This can be 
shown as follows. Let us consider the time development of a functional of e^^^, 
F[e^^^]. The time evolution equation is given by 



In this expression, the contribution of the auxiliary fields defines the infinitesimal 
local gauge transformation. The generator of the infinitesimal local gauge transfor- 
mation is defined by 



This satisfies the supersymmetric extension of the Gauss law algebra defined by 



Therefore, the introduction of the stochastic gauge fixing term has no effect on the 
time development of local gauge invariant observables, while it breaks the extended 
local gauge invariance, (3-10) and (3-11). This is the reason why the self-interaction 
of the vector superfield in the stochastic gauge fixing term simulates the Faddeev- 
Popov ghost contribution in the path-integral method. 

Here we comment on the implication of the extended local gauge transformation, 
(3-10) and (3-11). In the weak coupling region, (3-11) reads (f) (l)+(3~^ ij+ig[<p, S] . 
The variation indicates a covariant derivative of the gauge parameter U with respect 
to the stochastic time. In this sense, we may interpret the extended local gauge 
transformation as a 5-dimensional local gauge transformation for which the stochastic 
time is identified with the 5-th dimensional coordinate and the chiral and the anti- 
chiral superfields play the role of the 5-th component of the gauge field. For the 
YM case, this 5-dimensional extension of the local gauge transformation in SQM 
was discussed in Refs.lO), 14), 15) and 17). We note that the extended local gauge 
symmetry in (3-10) and (3-11) is not a 5-dimensional supersymmetric extension. The 
role of the stochastic time is distinct from that of the other 4-dimensional space-time 
coordinates, and the global supersymmetry is the 4-dimensional one. 





(3-15) 



G(0i, </.l),G(02, h) =-iG{[cj)i,cj)2], [</>!, </'2]) . 



(3-16) 



§4. Equivalence to the Faddeev-Popov prescription 



In the YM4 case, as we explain in §3, the equivalence of the stochastic gauge 
fixing procedure and the Faddeev-Popov prescription in the path-integral method is 
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shown on the basis of the BRST invariance of the Faddeev-Popov probabihty dis- 
tribution (3-5).^)'^) Here we give a proof of the assertion that the stochastic gauge 
fixing procedure proposed in the previous section for SYM4 in the superfield formal- 
ism reproduces the Faddeev-Popov distribution given in the path-integral method. 

To show the equivalence, we first define the Faddeev-Popov distribution of SYM4 
in the BRST invariant path-integral method. Let us introduce the BRST operator 

^{~LV^ + ^/L,')J^}, (4.1, 
where [X x y]" = fo-bcj^byc fpj^g Faddeev-Popov distribution is then given by 

Jgh 

'S'tot = S + Sgf + Spp , 

Sgi + Sfp = JcfzSBRSTF + TV{ Jd''zBf + jcffW] + 2^ j(f zj f] , (4-2) 

where J^^ = J VBVBVcDcDdVclVfVf. Here = {B,c,d,f) and (p = {B,c,cf,f) 
are chiral and anti-chiral superfields, i.e., they satisfy Dp = Dp = 0. / and / are 
variables for the gauge averaging. The function F for the gauge fixing is defined by 

F = -i j(fz{d + (^fV . (4-3) 

The partition function is given by -Z^ = JDVPyv in the BRST invariant path-integral 
formulation. 

We show that the Faddeev-Popov distribution, Pfp given by (4-2), is a stationary 
solution of the gauge fixed Fokker-Planck equation 

+^(^-e2^^^<^))V[r,e2«^] , (4-4) 

which is derived from the stochastic-time evolution equation of the expectation value 

of observables in (3-14) under the appropriate choice of the gauge fixing functions 
(f) and (f). Namely, the r.h.s. of the Fokker-Planck equation (4-4) vanishes for PpPj 
provided that the gauge fixing functions are defined by 



r{z) ^ Ag^ J ^Jd'z'S'{z){si''\z'){£i''\z')F)e-'^°^]p^^ . (4-5) 
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The proof is as follows. Substituting the distribution Ppp into the Fokker-Planck 
equation (4-4), the r.h.s. reads 



+^^{4>{z)-e'9^y<P{z)rPpp} . (4-6) 

The first term in (4-6) reads 

Jd'zSi^\z){gr'^^c^{£^^\z)F)+6BKSTi£i''Hz)^^^ 

-\ fd'^Sf\^){9r"'o\^f\^)F) + W(^i^)(^)F)}e-^- , (4-7) 
while the second term in (4-6) is evaluated as 

ig Jd^zd^z'^{z)£i^^ {z') {Si^^ {z')F)e-^'°' } 

-ig j^^ jd^zd^z'c-{z)£f^ {z){Ef^ {z'){8f^ (/)i^)e-^*°* } 
= \j^^ jd'z5BKST{£i''^ (z) (4^) (^)F)e-^-* 

= 11^^ jd'z[ [feRST, {z)] (^f ) (^)F)e-^- 

+ [feRST, Sf^ {z)] {Ef^ (^)F)e-^- } 

+\ j^^ ld'z{si^\z)5BKST{ii'^^ (z)F)e-^- 

+^(«) (^)5brst (4^) (z)F)e-^*- } 
= - the first term in (4-6) . Q.E.D. (4-8) 

In the above proof, we have used the commutation relations 

[W, ii'^\z)] = -gr"^c*{z)i^,'^Hz) , 

[5brst, ii^'Hz)] = -gr'^c\z)if\z) . (4-9) 

Hence, the distribution Ppp is a stationary solution of the Fokkcr-Planck equa- 
tion. However, as we have already mentioned, this does not imply the existence 
of the unique equilibrium limit: limT-^^ P{t) = Ppp. To demonstrate the relax- 
ation of the probability distribution to a unique equilibrium distribution, we need 
a regularization procedure in order to reduce the number of the degrees of freedom 
that preserves the global supersymmetry as well as the local gauge symmetry. Fur- 
thermore, if the Gribov problem exists in the covariant gauge fixing (4-3) in the 
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conventional approach, the positivity of Ppp is broken. Then, it is not justified that 
the probabihty distribution at finite stochastic time P{t) uniquely relaxes to the 
Faddeev-Popov distribution Ppp in equilibrium. We leave the investigation of these 
delicate problems to future works. 

The gauge fixing functions (4-5), which are essential for the proof of the equiva- 
lence, are different from those in (3-12), which are convenient for actual calculations. 
The proof of the equivalence relies on the particular form of the gauge fixing functions 
given in (4-5). Before ending this section, we comment on the relation between these 
two gauge fixings, (4-5) and (3-12). The gauge fixing functions (4-5) are expressed 
as 

cl>^ = 2ig Jd'z'{c%z){c' + c')\z'))L%{z')£i''\z'){S + S'^,) 
-2ig j<fz'L\{z')£f\z'){c%z){c' + c^)\z')) , 
= 2ig Jd'z'{cf{z){c' + ci)\z'))L,\z')Si'^\z'){S + S'^) 

-2ig jd^z'L,\z')iP{z'){<f{z){<^ + f)\z')) , (4-10) 

where S'^ = Jd^ z- {D'^V)"' {D'^V)"' , which is obtained by the integration of the aux- 

8 _ _ 

iliary superfields B, B, f and / in the definition of Ppp (4-2). The ghost superprop- 
agators in (4-10) are defined by Ppp, for example, 

{Ca{z)c'\z')) ^ Pp-p^ / Caizy\z')e-^*''^ . (4-11) 

Jgh 

We can expand each quantity in (4-10) with respect to the coupling constant g ( in 
other words, with respect to the vector superfield V ). In principle, it is possible 
to evaluate (4-10) pcrturbatively with the appropriate renormalization procedure of 
the singularities. In order to determine the leading term, we substitute the free 
propagators of the Faddeev-Popov ghost, 

= 2{-\D\{-\D\,{a,z\{-U)-^\z',h) , (4-12) 

into (4-10). Then, we obtain 

<i) = i^D'^D'^V + 0{V'^) , 

^ = -i^D'^D^V + 0{V'^) . (4-13) 

This is what we alluded to in the previous section. Namely, the local stochastic gauge 
fixing functions (3-12), which are almost unique from the dimensional analysis, are 
the linear terms in the stochastic gauge fixing functions that are equivalent to the 
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conventional Faddeev-Popov prescription. This fact also implies that the stochastic 
gauge fixing term with (3-12) introduces a drift force for the longitudinal mode of 
the gauge field in the superfield Langevin equation (3-9). In general, as is clear from 
(4-10), the local gauge fixing in the conventional approach yields non-local stochastic 
gauge fixing functions. On the other hand, if we chose the stochastic gauge fixing 
functions to be local, the corresponding gauge fixing would be non-local in the 
conventional approach. As we have shown in (3T4), the expectation values and 
the correlations for local gauge invariant observables are independent of the gauge 
fixing functions, and the difference between these two gauge fixings is irrelevant in 
a perturbative sense. 

§5. Stochastic action for SYM4 and the BRST symmetry 

For perturbative analysis in SQM including the renormalization procedure, it 
is convenient to introduce the so-called stochastic action by the path-integral rep- 
resentation of the SQM approach. The construction of the stochastic action is 
standard. For SYM4, applying a general formulation for the first-class constraint 
systems,^^^ the BRST invariant stochastic action is defined as follows. 

We first consider the continuum limit of the Langevin equation (3-9) and the 
noise correlation (2-1) by taking At 0. This limiting procedure defines the con- 
tinuum Langevin equation E{V{t, z)) = r]{T, z), where the noise superfield is defined 
by ri{T,z) = liniAr^o -^"^iT, z)/ At, and its correlation is given by the continuum 
limit of (2-1). E{V) is defined by the continuum limit of (3-9). Then we insert unity 
in the form 

1 = I VGVVV(I)V4>6{E{V) - r/)det(^^)(5(#)(5(#) (5-1) 

into the Gaussian integral representation of the correlation of the noise superfield ry. 
Here, the functions # and # are defined by 

^ = - ^^L>2^V , ^ = (j)+^(,D'^D'^V , (5-2) 

for the stochastic gauge fixing (3T2). In the Ito stochastic calculus, we choose 
6{0) = 0, so that the dynamical variables arc not correlated with the equal stochastic- 
time noise variables for the solution of the Langevin equation. This gives a trivial 
Jacobian, det{6E/6V) = 1, while a careful evaluation of the integration measure^^^ 
yields the contribution 

exp(l/d»*G-(.,.)j^^) , (5.3) 

which includes (5^(0) in this case. Employing the regularization by dimensional re- 
duction,'^^)' we discard this contribution. 

The integral representation of the (^-functional in (5-1) requires an auxiliary 
superfield tt" that is not a vector superfield. By integrating out the noise superfield 
ri, we obtain a stochastic action as a functional of the vector superfield F", the 
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auxiliary superfield tt", and the chiral and anti-chiral superfields cp and (p. After 
the redefinition of the auxihary field tt" = L'^i^w^, w°' becomes a vector superfield, 

As we have shown in (3-10), this stochastic action is also invariant under the 
extended local gauge transformation, provided that the auxiliary superfield vr is 
transformed as tt" e^'^'^^t-Tf". This also induces the local gauge transformation 
of the auxiliary vector superfield w. Therefore, we can introduce a BRST trans- 
formation by making the replacement X! —Ac, E'^ —Ac. Here, A is a purely 
imaginary Grassmann number. The BRST symmetry is realized on the field con- 
tents {(F", Wa), (0",S",c»,Ca), (0",^'',c'^,a„)}. Here, (/)'^, 5", c« and P are chiral 

— - a 

and (f)"', B , and Ca are anti-chiral superfield. As is clear from (5T), we impose 
the gauge conditions (3T2) with the Nakanishi-Lautrup fields, B and B, and the 
functions, ^ and ^. This choice of the gauge conditions is an extension of the "flow 
gauges "in the YM case, in which the Faddeev-Popov ghost decouples from the 
gauge field. Explicitly, the extended (or 5-dimensional) BRST transformation 
is defined by 





= -A^nc" + A^c'^V , 






Sbkst(P"' 


= A2kc" + Xg[(t) X c]" , 


f^BRST^" 


= A2kc° + Xg[4> X c]" , 


<5brstc" 


= -A|[cxc]«, 


i^brstc" 


= -A|[cxcT, 


<^BRSTCa 


= iXBa , ^BRSTCa = iXBa 


(^BRST-B** 


= , (^BRST^" = . 



(5-4) 

Here, we have redefined the scaling parameter of the stochastic time as (3 = — — . 

2k 

(Also see Appendix A for the convention.) The extended (or 5-dimensional) BRST 
invariant stochastic action is given by 



(^^e^BRST ^ = VVVwVcDcDcDcD(t)V(f)VBVB , 

K =Jd^zdT ^j-G-b^aWb + iWa{v'' - ^ (^(P'^W^)' + (^dW^")'V) } , 



4k 



*^ The hermiticity of Wa is defined only from the consistency of the transformation prop- 
erty. In analogy to the Ostarwalder-Schrader conjugation in the Euclidean space-time, the pre- 
cise definition of the hcrmitian conjugation includes the time reversal with respect to the stochastic 
time,^^^ Viry'' = V^-t)", w(r)l = zuai-T)-idS/6V''{-T), for the hermiticity of the stochastic 
hamiltonian, H = i fi,^ zrtUaV'^ — K. 
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-^BRST = K + K]SIL + Kpp , 

i^NL + Kfp = -i- j(fzdT[s(ca^'' + 5a^«) } . (5-5) 

Here, ^brstO = \{sO). <P and # are defined by (5-2). This formulation is sufficient 
to perform the perturbative analysis. As we explain in the following, we may consider 
truncated versions of (5-4) and (5-5). 

We redefine the Nakanishi-Lautrup fields, B and B, and the anti-FP ghosts, c 
and c, as 

b=(--dAb. b=(--dAb 



<^^{-\d')1, -d^(-\Dy. (5-6) 
The BRST transformation for anti-FP ghosts becomes 

(^BRSTc'a = i\Ba , feRSTcJj = iXBa ■ (5-7) 

Correspondingly, the BRST exact part, -RTnl and i^pp in (5-5), is expressed in terms 
of the chiral measure in the superspace, d^z = d'^xcfd and d^z = d^xd^O, as 

i^NL + = {K j<^^dTC^a^^ + jd!^ zdTc'^0''') } . (5-8) 

To reduce the number of auxiliary fields, we consider a consistent truncation 
of the BRST transformation (5-4) by integrating out the auxiliary superfields, ^ 
and (j), and the Nakanishi-Lautrup superfields, B and B. Since cp and (p are chiral 
and anti-chiral, respectively, the integrations of these superfields yield 5-functionals 
which impose the conditions 



iBa = {-\d^) (l^bLa") - 9[C X c']a , 



iBa = -(-]d^) (l^bL\) - g[c X c']a . (5-9) 



4 J\2 

The equations of motion (5-9) yield a truncated BRST symmetry, as well as 
a corresponding truncated BRST invariant stochastic action. The truncated BRST 

transformation is expressed in the extended phase space | (F", Wa), (c", c"*), (c", c'")| 



as 



5brstF« = -A|lV' + a|c*V , 

feRST^a = X-j^daL^b'^cC^ ~ ^-j^daL^^Wc , 
^BRSTC" = -A| [cxcf , 

5brstc" = -a|[cxc]" , 

^BRSTc'a = -X(-^D'^^ {^WbL^^ - \g[c X c']„ , 

Wc4 = {\^a^b) - M[C X Aa ■ (5-10) 



'short' title 



15 



We note that the truncated BRST transformation is nilpotent. It preserves the 
chirahty of the superfields; i.e., Dc = Dc/ = Dc = Dc' = 0. The truncated BRST 
invariant stochastic action is reduced from the 5-dimensional extended one (5-5): 

Z = j dp.e^^^'^ , dii = VVVwVcDc'VcDc' , 

Kbrst = K' + Jd^zdTc'a(f + Jd^^zdrc'^c" 

fd^zdrlsic'^D'^V'' - a'aD^V")] . (5-11) 
2k J I J 

Here, ^brstC = X{sO), and K' is given by K with = = in (5-5). The 
BRST invariance of this truncated stochastic action is manifest, if we regard (5-11) 

as a Legendre transformation, because of the invariance of the derivative terms with 
respect to the stochastic time. Namely, we have a non-trivial BRST cohomology of 
(5-10); (5BRST|i fd^zdTWaV'' + jd^ zdrc' ac"- + Jd^zdrc'^t^ = 0. We also note that the 
canonical conjugate momentum of the FP ghost c (c) with respect to the stochastic 
time is c' {(f). The BRST exact part in (5-11) is decomposed as 



^^ jd^zdrd'^^^^icf.D^V^-c'aD^V-) 



jd^zdrv^ai^KD'^D'^V^ + L%D^D^V''^ 

jd^zdT[\(sf,D^ - c',D'){L\c' - c^V) 

-g[c X cf]aDW + g[c x c']„L>V"} . (5-12) 

The first term on the r.h.s. of this expression is the stochastic gauge fixing term. The 
second and third terms describe the interactions between the vector superfield and 
the FP ghosts. We note that the truncated BRST symmetry (5-10) is independent 
of the specific form of the stochastic gauge fixing functions. For example, it allows 
the gauge fixing functions given in (3-13). 

The stochastic action (5-11) is also invariant under the scale transformation of 
Faddeev-Popov ghosts: 

c ^ e^c, c' e-Pc', c e^c, c' e-^c' . (5-13) 

We assign the ghost number +1 for c and c, — 1 for c' and c', and for the remaining 
fields. 

In order to define the superpropagators, we determine the kinetic term by taking 
the limit 5 — > 0. It is given by 

^BEST =J<^^<'r bk^''^' + iro„{l>° + i ( - □ + i^(D'D^ + C^D^)) V}] , 
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For example, in the Feynman gauge, ^ = 1, we obtain the simplest superpropagators: 

^ ^ ^ V ' // k{ilu + k^/K){-tUJ + k^/K) ^ ^ ^ ' 

{x^\t,z)w\t',z')) = Q, 

(c-(r, z)-c\r', z')) = - ^.^ " ^0^^(^" - ■ (5-15) 

Here, we have suppressed the momentum integration Jd'^kdui{2Tr)~^e^'''^^~^'^'^^'^^'^'''^'\ 
We have also defined the covariant derivatives Dq, = da — c^A^"^"! ^^'^ = 
—da + G'^(^^akm- As is clear from these expressions in (5-15), the superpropagators 
for the vector superfields, V and w, are the same as those for the YM case, except 
the (5-functions for the Grassmannian coordinates, while the superpropagators for 
F-P ghosts include the chiral projection operators, and D^. We comment on the 
retarded (or causal) property of the superpropagators for (Vzu) and the F-P ghosts. 
In the superfield Langevin equation, we have to choose 6{0) = 0, so that the dynam- 
ical fields are not correlated with the equal stochastic-time noise fields in terms of 
the Ito stochastic calculus. This prescription is also a prerequisite in the definitions 
of the stochastic action and the superpropagators. Therefore, in the configuration 
space, the choice ^(0) = is always assumed for the definition of {Vzu) and the ghost 
sup erpropagator s . 

The superpropagators and the BRST invariant stochastic action (5-11) define 
the perturbativc expansion with supcrgraphs in the SQM approach by preserving the 
global supersymmetry in the superfield formalism. We note that we have specified 
the regularization procedure by employing regularization by dimensional reduction 
to discard the singularities, such as S^{0). 



§6. Power counting and the renormalizability of SYM4 
in the SQM approach 



In the path-integral approach, the perturbativc analysis is carried out in terms 
of the supergraphs. This ensures the non-renormalization theorem in SYM4.^^)'^^^ 
As we have explained, the perturbative calculation with the supergraphs can also 
be applied by using the superpropagators defined in (5-15) in the context of SQM. 
Therefore, the non-renormalization theorem also holds in the SQM approach, which 
is a consequence of the fact that the non- vanishing contributions in the perturbative 
expansion come from the loops which include the four spinor covariant derivatives 
DaDfiDaDp to eliminate 5'^{e - 9')d^{0 - I (?=(?', 9=e' in the superspace 5-function. 
Furthermore, the symmetries and the simple power counting argument indicate the 
multiplicative renormalizability of the SQM approach in terms of the BRST invariant 
stochastic action for SYM4. 
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Let us estimate the divergence of a Feynman diagram G{V, I, E) with V vertices, 
/ internal hnes, E external lines and L = I — V +1 loops. More precisely, we define, 

E^ + Ey ■ 

# external lines of V 

# internal lines of (V^V^) 

# vertices with two w and nV.{Q-l) 

We also have the topological relations 

^(^1." + ■^s.n) = + Iv^ + Ev, Y.^Vi,n + 2V2,n) = Iv^ + E^. (6-2) 

n=2 n=2 

The Feynman diagram G(y, 7, E) may be evaluated as 

G{VJ,E) 

Thus the degree of the ultra-violet divergence of the Feynman diagram G{y,I,E) is 
estimated as 

4L - 2Ivv - 2/ + 2 Vi,n) =4.-2E^. (6-4) 

n=2 

It is necessary to introduce two types of counter-terms to cancel the divergences, 
which are constrained by the Ward-Takahashi identities of the BRST symmetry 
(5-10): 

(i) For E^ = 2, the logarithmic divergences are cancelled by terms of the form w'^V'^, 
with n = 0,1, 2, ...oo. 

(ii) For E^ = 1, the logarithmic divergences are cancelled by terms of the forms, 
such as wV and zuD°'D^Dc(V'\ with n = 1, 2, ...oo. 

(iii) For E^ = 0, there are no Feynman diagrams. 

We must pay a special attention to the E^ = case. Since the superpropagator 
{Vzu) is a retarded one with respect to the stochastic time, the contribution of a loop 
vanishes if it consists of only the superpropagator (Vw) as the internal lines. This 
means that the Feynman diagram includes at least one external line of w. There is 
no ambiguity due to the regularization procedure with regard to this fact, because we 
have to choose ^(0) = in the Ito stochastic calculus. This is an essential point for 
the renormalizability of the present approach. Without this property, it is obvious 
that the SQM approach is not renormalizable. The same conclusion is obtained from 
the stochastic Ward-Takahashi identity. -"^^^ We note that, for E^ = 1, the divergence 



I = Ivvu + Ivv, E = 

E^ : # external lines of zu, Ey '■ 

Ivvj '■ # internal lines of (yw), lyv '■ 

Vi,n '■ # vertices with w and n V, F2 n : 
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appears to be quadratic. However, the BRST invariance reduces the order of the 
divergence to logarithmic. 

The power counting of the Feynman diagrams with FP ghost external lines is 
similar to what we have explained for the Feynman diagrams with only external 
vector superfields, except that the ghosts are the chiral and anti-chiral superfields. 
The degree of the divergence is given by 

4-2E^-^{E^ + E^,+E^ + E^,), (6-5) 

where Ef., E^', E^ and E^' are the numbers of the external ghost lines for c, c', c and 
c', respectively. We have also used the ghost number conservation to derive (6-5). 
Because of the retarded (or causal) nature of the FP ghost superpropagators in 
(5-15), the ghost loops do not contribute to the renormalization of the 1-P-I vertices. 
This means that, in the actual calculation of physical quantities such as the [3- 
f unction of the gauge coupling.^ we m,ay discard the F-P ghost sector. On the 
other hand, since the vector superfields couple with the ghosts and they contribute 
to the renormalization of FP ghost sector, (5-12) and Jdr^Jcfzc! ac"" + Jd^zclJ}^) 
are renormalizcd by the loop correction of the vector superfields. Therefore, the 
existence of the BRST symmetry is necessary for the formal proof of the manifest 
renormalizability of SYM4 in this context, in order to apply the standard arguments 
in terms of the Ward-Takahashi identity. As we have already noted with regard 
to (5-11), the non-trivial BRST cohomology (5BRST{/(i^r(i jd^zWaV" + jd^zc'ac" + 
Jd^zc'g^cf)} = plays a key role for the renormalizability of the F-P ghost sector, and 
this is the reason why we have introduced the scaling parameter k. The multiplicative 
renomalizability requires this additional coupling constant k, due to the existence of 
this particular non-trivial BRST cohomology. Since the stochastic gauge fixing term 
and the remaining F-P ghost term takes a BRST exact form, the power counting 
argument, the BRST symmetry, and the ghost number conservation indicate that 
the stochastic action (5-11) is multiplicatively renormalizable by the wave function 
renormalizations and the renormalizations of the couplings g, k and ^. 

We add a comment concerning the proof given in Appendix B of the renormal- 
izability of the SQM approach in terms of the BRST invariant stochastic action. In 
the truncated version (5-11), since the truncated BRST symmetry is independent of 
the stochastic gauge fixing term, there is an ambiguity for the renormalization of the 
gauge fixing term, which cannot be determined by the symmetry argument. As we 
have already explained as an example, we can change the stochastic gauge fixing, 
replacing V by h{V) with an arbitrary function h, without breaking the truncated 
BRST symmetry. Therefore, in a formal sense, this ambiguity of the renormaliza- 
tion for the gauge fixing term could break the multiplicative renormalizability of the 
gauge fixing and the F-P ghost sector in the truncated version (5-11). In Appendix 
B, we prove the renormalizability of the BRST invariant stochastic action principle 
in terms of the Ward-Takahashi identities for the 1-P-I vertices. To avoid the ambi- 
guity in the renormalization of the gauge fixing term, we prove the renormalizability 
of the extended BRST invariant stochastic action (5-5) instead of (5T1). We obtain 
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the renormalized stochastic action and the counter-terms as 

1/2 

Co) = i^BRST( ^99, Z^K, Z^i) , (6-6) 

where = {V,w,c,c,cf,(f,(l),^,B,B). The subscript "0 "denotes the bare quanti- 
ties. The renormahzation constants satisfy the relations 

^-14/'^^' = 1, = ^^4/' = Zb = Z^ = Z,, = Z,, . (6-7) 

The first two relations in (6-7) arc consequences of the non-renormalization of the 
stochastic gauge fixing term in (5-5). 

Since the truncated version is derived from the extended version by integrating 
out the auxiliary fields {(f>,(l), B, B), the two formulations are formally equivalent. 
Furthermore, the expectation values of local gauge invariant observables are inde- 
pendent of the stochastic gauge fixing functions, and F-P ghosts do not contribute 
to the renormahzation. Therefore, the ambiguity in the renormahzation of the gauge 
fixing term in the truncated BRST invariant formulation may be harmless for actual 
calculations. In fact, the perturbative analysis in terms of the truncated version 
of the stochastic action in the background field method is consistent with that of 
the conventional approach in the one- loop approximation, which is confirmed by the 
explicit calculation of the /^-function of the gauge coupling.^^^ 

§7. Discussion 

We have introduced the stochastic gauge fixing procedure for = 1 SYM4 
in the SQM approach in which the superfield formalism preserves the local gauge 
symmetry as well as the global supcrsymmetry in the underlying stochastic pro- 
cess in the sense of Ito calculus. Wc have shown that the standard Faddccv-Popov 
probability distribution is a stationary solution of the Fokker-Planck equation, and 
it is conjectured to be reproduced in the equilibrium limit. This means that the 
contribution of the Faddecv-Popov ghost superfield in the path-integral approach 
is simulated by the extra self-interaction of the vector superfield in the stochastic 
gauge fixing term. The perturbative properties have been studied in terms of the 
BRST invariant stochastic action with regularization by dimensional reduction. In 
particular, the non-renormalization theorem holds, as in the conventional approach. 
Wc have also shown the perturbative renormalizability of the present formulation 
based on the power counting argument and the BRST symmetry of the stochastic 
action. It is essential for the SQM approach to preserve the manifest supcrsymme- 
try of SYM4 that the supcrsymmetry is defined as a translation in the superspace in 
terms of the superfield formalism. In the SQM approach, a possible breakdown of 
the supcrsymmetry would be induced by the correlation of the noise variables. In the 
superfield formalism, the correlation of the noise superfield is also invariant under 
the translation in the superspace. Therefore, the time development with respect to 
the stochastic time preserves the supcrsymmetry, and thus the system is invariant 
under the translation at finite stochastic time. In this respect, SQM can be ap- 
plied directly to the supersymmetric models formulated in terms of the superfield 
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formalism. By contrast, if we apply SQM to a supersymmetric model that consists 
of only physical component fields, we would encounter an artificial breakdown of the 
supersymmetry at finite stochastic time, though the supersymmetry is expected to 
be recovered in equilibrium. 

In this paper, we have formulated the stochastic action on the basis of the BRST 
symmetry. Since the interaction is non-polynomial for SYM4 in the superfield for- 
malism, it is convenient to apply the background field method (BFM) for pcrturba- 
tive calculations in the SQM approach. In the conventional (path-integral) BFM 
for SYM4 without chiral matter multiplets,^^) the chiral superfields, such as the 
Nakanishi-Lautrup fields and Faddeev-Popov ghosts, are defined by the chiral condi- 
tion with respect to the background covariant derivative. This causes a non-trivial 
contribution of the Nakanishi-Lautrup fields. In order to cancel this contribution, 
we need an additional ghost, i.e., Niclscn-Kallosh-type ghost. In Ref.37) we car- 
ried out a one-loop calculation to determine the /3-function of the gauge coupling by 
applying BFM to the stochastic action. The one-loop /3-function is identical to the 
known results in the path-integral approach. This indicates that the contribution of 
the Nielsen-Kallosh ghost can be simulated by defining the stochastic gauge fixing 
function in terms of the backgound covariant derivatives, instead of and in 
this paper. We note that the necessity of the Nielsen-Kallosh ghost is peculiar to 
the conventional BFM. 

One important check for the equivalence of the SQM approach to the stan- 
dard path-integral approach is to evaluate anomalies in SYM4.^*) In Ref.39), we 
derived the anomalous Ward-Takahashi identity for the superconformal symmetry 
in this context with the stochastic gauge fixing procedure in the one-loop approx- 
imation. We also confirmed the one-loop equivalence of the effective equation of 
motion obtained in the SQM approach and the conventional one, which includes the 
contribution from ghosts. The equivalence in BFM is shown by expressing the ghost 
superpropagators in terms of the vector superpropagator with the help of Slavnov- 
Taylor identities and Schwinger-Dyson equations at the one-loop level. This explicit 
one-loop equivalence in BFM also supports the formal proof given in this paper. 
In the conventional approach, there is no anomaly for the local gauge symmetry, 
because the gauge multiplet and the ghost multiplet belong to the adjoint represen- 
tation. This is also true in the SQM approach, due to the one-loop equivalence in 
BFM to the conventional approach. In this paper, we have considered only the vector 
superfield, i.e., the gauge multiplet in SYM4 including F-P ghosts in the conventional 
sense. The extension to introduce matter fields, such as chiral matter multiplets in 
any representations, does not appear to present any obstacles for the case of the 
anomaly free combination. 

We have mainly investigated perturbative properties of SYM4 in the SQM ap- 
proach with the stochastic gauge fixing procedure. In order to give a precise definition 
of the stochastic process for SYM, it is necessary to reduce the number of degrees 
of freedom to a finite value for the existence of a unique equilibrium limit. The best 
way to carry this out may be to apply the lattice regularization. However, in the 
lattice regularization, the exact global supersymmetry is not manifest, though it is 
expected to be recovered in the continuum limit. In this respect, we may regard 
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the naive zero volume limit of SYM at large A'' as a regularization procedure that 
preserves the global supersymmetry and allows us to define the stochastic process 
precisely, because the finiteness of the partition function and correlation functions in 
this limit has been proven with finite N.'^^'^ In addition, though the fcrniion integral 
in SYM yields a Pfaffian, it is positive semi-definite in the four-dimensional case. 
Therefore, the zero volume limit of SYM4 in the superfield formalism, which is ex- 
pected to recover the planar limit of the continuum SYM4, may provide a toy model 
to study the stochastic process with exact global supersymmetry in the context of 
the SQM approach. 

We finally comment on some difficulties in the application of SQM to the con- 
structive definition of superstring theories in terms of supersymmetric large N re- 
duced models. As mentioned above, without auxiliary fields in the superfield 
formalism, we would encounter an artificial breakdown of the supersymmetry at 
finite stochastic time. Although we are able to expect that the supersymmetry is 
recovered in the limit of infinite stochastic time for the cases of SYMio and its zero 
volume limit, there may be no merit of considering a naive application of SQM to 
such models, because the Fokker-Planck hamiltonian, which is Lorentz invariant by 
construction, does not commute with the generator of the supersymmetry. One way 
to avoid this problem is to integrate out the fermion fields. However, in contrast to 
the case for SYM4, the Pfaffian is complex in the ten-dimensional case, and this may 
cause another difficulty in attempting to define the stochastic process in terms of 
the Langevin equation. These problems remain to be solved, and we hope that the 
present analysis of SYM4 in the superfield formalism provides an insight for future 
application of SQM to supersymmetric models, such as SYM in the Wess-Zumino 
gauge, which consist of only physical component fields. 
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Appendix A 

Conventions for SYM4 in the SQM Approach 

In this appendix, we summarize the superfield formalism of SYM4 in the SQM 
approach. We use a convention closely related to that of Rcf.30). In particular, 
the reality condition on the vector superfield V"^ = V is defined in the Minkowski 
space-time with the signature (— , -|-, -|-, -|-). In order to assure the convergence of 
the superfield Langevin equation, we consider the Wick rotation to the Euclidean 
space-time by ix^ = x^. In this case, the generators of the SL{2, C) are mapped to 
those for SO{4) = SU{2) x SU{2) by preserving the relation a'^dm = (^(^)9m^- It 
is not necessary to change the algebraic rule for the D-algebra, which is obtained 
by replacing ry™" grnn_ example, □ = 9f + c?! + c?! + c^l is understood. In 
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the Euclidean space-time, however, since 9a ^ Oa, the superfield formulation needs a 
doubled superspace, {x,9a,0d) = >S'i and {x^,9l^,da) = S2, where 9a and are two 
independent spinor coordinates.^^^ Thus, the number of degrees of freedom of the 
vector multiplct arc also doubled, and hence Aq 7^ A^, Vm 7^ v"^, etc. In other words, 
suppose that the vector superfield is defined on ^i, i.e., V G ^i, the hermitian 
conjugation is defined on ^2, i.e., G iS'2. Then the hermiticity of the vector 
sTipcrficld is defined by the self-conjugation in the sense of Ostewarlder-Schrader by 
identifying the doubled superspace, Si and 5*2, with the simultaneous time reversal 
operation. The identification of the doubled fermionic component fields is defined by 
the Majorana condition, Aa = A^, by continuing back to the Minkowski space-time. 
This prescription is also applied to the chiral superfield. If we difine a chiral superfield 
and an anti-chiral superfield 'F on Si, i.e., D^'I^ = Da"^ = 0, the hermitian 
conjugation is defined on 5*2, i.e., •f'^ G ^2. Then can be identified with # G 5i in 
the sense of the 0-S conjugation. Therefore, if we study the hermiticity of the vector 
superfield in the Euclidean superspace in terms of this doubled formulation, we obtain 
no physically distinct conclusions in comparison with the hermiticity condition of the 
vector superfield in the Minkowski space-time. In this respect, it is safe to use the 
superfield formulation of the vector superfield in the Minkowski space-time, even 
after the Wick rotation, if wc keep in mind that the hermitian conjugation is defined 
in the sense of Ostewarlder-Schrader. 

The vector superfield V is SU{N) algebra valued: V G su{N) ; V = VHa, 
[ta, tb] = ifabctc and Tr{tatb) = kSab- For the convenience, is assumed to be in the 
fundamental representaion. The action for SYM4 is defined by 




Here 

Wa = - \D^e-^^^Dae^^^ , Wa = lo'e'^^Dae-'^^ . (A-2) 

The analogues of the Lie derivatives and the Haar measure defined on a group 
manifold are introduced as follows. The differential operator si^' (z) is defined by 



Here 



25^" ^^hvb{z) ' 
2gL^ 



(A-4) 
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where z = ix,9,9), 6^{z-z') = 6^{x - x')6^{9 - e')S'^{e - 9'), LyX = [V,X],^^^ and 
= La K These differential operators satisfy the algebra 

[£i^\z),£j,'^\z')]=-if,,,S'{z-z')S(^\z) , 

'ii''\z),£l''\z')] = ifabcS'iz - z')i(/\z) , (A-5) 

which imply 

L'^a^cL'^b ~ ^'^b^cL'^a = '^'^9 fabc^'^c i 

L:dcL\-L%dcLj = 0. (A-6) 

We also define a measure by introducing the metric = L^{z)L^{z) = 

L\{z)L\{z) and Gab{z) ^ K,-{z)K,-{z), where K,%z)L,^{z) = L,%z)K,\z) = S.K 
By the metric, the integration measure is defined by y/GW . 
The superfield Langevin equation for SYM4 is given by 

-^(Z\e25^)e-25^(r, z) = -pAr2g£^^\T, z)S + Aw{t, z) . (A-7) 

Here ^W(r,z) = tj^^^ {t,z). The equations of motion are defined by VW = 

Wa}, VW = ^Va, VF"}, where = e-^^V jj^^2gV ^^^^ -p^ = ^2gV ]j^^^2gV ^ 

Then (A-7) becomes (2T). Since the reality condition implies e^^^^^PW = PW,^^) 
we may discard one of them in (2-1). 

The probability distribution P(r, e^^^) is defined in terms of the expectation 
value of an arbitrary observable F[e^3^] as 

(F[e2»^](r)) = jF[e^sV^PiT,e'^3V^^jyy _ (^-8) 

Here, the expectation value on the l.h.s. is evaluated, using the solution of the 
Langevin equation (21), by accounting for the correlations of all the noise superfields 
{Auj\Alo{t') < t - At}. 

The Langevin equation (A-7) (i.e., (2-1)) is converted to that for the vector 
superfield V which appears in a geometrical expression 

AV^ir, z) = -PAtG''\t, ^)^^ + A^S\t, z) , 

{A^E\t, z)A^E\t, z'))^^^ = 2(3AtG'^\t, z)d'{z - z') . (A-9) 

Here, AyjE"- is a collective noise superfield defined as AyjE"'{T, z) = Aw^{t, z)L^{t, z). 
In the matrix representation, (A-9) becomes (2-4). In terms of the Ito calculus, the 
noise variable is uncorrelated with the equal stochastic-time dynamical variables. 
This means that the r.h.s. of the collective noise correlation (A-9) is not the expec- 
tation value. 
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The vector superfield must remain real under the time evolution described by 
(A-7), i.e., V{t + Ar)"^ = V(t + Ar). To satisfy this constraint, we impose the 
hermitian conjugation of the noise superfield as 

Aw{t)^ = e-'^3^^Aw{T) . (A- 10) 

This is consistent with the transformation property of Aw given in (2-3). Although 
Aw is not a vector superfield, it follows from (A- 10), 

A^E-^ = Aw'^L,''^ , 

= '-Tr( Aw^ . t^) = A^E^ . (A-ll) 

Thus, the collective noise superfield is a vector superfield. This ensures that the time 
development described by (A-9) (i.e., (2-4)) preserves the hermiticity of the vector 
superfield. 

The Fokker-Planck equation corresponding to (A-9) is expressed as 

^^(-.V')=;3/.».^{g-(.,.)(^ + ^)p(.,7)} . (A.12, 

The hermitian conjugation of the noise components is uniquely determined from 
the hermiticity assignment (ATO). Here we summarize the results in the "almost 
Wess-Zumino gauge " in which we can expand 

Aw^ = Aw - 2g[V\^^z, Aw] + 2g^[V\^z, [V\^z, Aw]] , 

y|wz = -Oa'^evmix) + io'^eXx) - ie^ex{x) + ^e^e^D{x) . (a-i3) 

This relation determines the imaginary part of the complex noise Aw. For example, 
for the vector component of Aw, (A-13) reads Aw\^^^^g = Aw\g^mg + 2g[vm, Aw\ ]. 
Here Aw\ is the lowest component of Aw, and Aw\^ = A^v\. This implies that the 
imaginary part of the vector component Awlg^mg is given by —g[vm, Aw\ ]. By 
contrast, the collective noise superfield AyjH has a manifestly hermitian expression, 

Ay,S = -Aw + -Aw^ + ^[F|wz, [l^lwz, Aw + Aw^]] . (A-14) 

Since AyjS is a vector superfield, it can be expanded into its components in a 
standard manner: AyjS = + ... — 6a'^6{AwS\g„rng) + ... . The correlation 

of each component of the collective noise is uniquely determined from (2-1). The 
details of the component expansion are given in Ref.l). For example, for the vector 
components, the correlation in (A-9) reads 

{{A^S\g^,no){T,x)iA^S\g^ng){T,y)) = -AfiAT5rnnS\x ~ x') . (A-15) 

In the superfield formalism, the precise meaning of a superfield expression is defined 
by its component expressions. In this respect, though the equilibrium distribution 
is independent of the value of the parameter (3, we need to choose /3 = — ^ in the 
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superfield Langevin equation and the noise correlation in (2-1) to provide a proper 
normalization in their component expressions with the conventions of Ref.30) in 
the Euclidean space-time. In this note, we introduce a scaling parameter k for the 

renormalization procedure. Specifically, we define /? = , as we define the BRST 

2k 

symmetry (5-4) and the BRST invariant stochastic action (5-5). 



In this appendix, we demonstrate the renormalizability of the SQM approach in 
terms of the BRST invariant stochastic action principle formulated in this paper. For 
the YM4 case, the renormalizability was proved formally in Rcf.l3). That proof is 
based on the Ward-Takahashi identities derived from the truncated BRST symmetry, 
which corresponds to (5-10) for the SYM4 case. The truncated BRST invariant 
formulation may be convenient, because of the reduced number of the auxiliary 
fields. However, the local identities for the changes of the integration variables 
c' c' + 6c' and c' c' + 6c' include composite operators that are not identified with 
the variation under the truncated BRST transformation (5-10). This fact complicates 
the application of the standard procedure to determine the divergent terms for the 
stochastic gauge fixing and the FP ghost sector in terms of the Ward-Takahashi 
identity. 

Instead of employing the truncated BRST formulation, we demonstrate the 
renormalizability in the SQM approach in terms of the extended (or "5-dimensional 
") BRST symmetry (5-4) for the extended BRST invariant stochastic action (5-5). 
Although, this includes extra auxiliary superfields ((/>, (j)) as the "5-th component 
"of the gauge field and the Nakanishi-Lautrup fields, {B, B), provided that we do 
not evaluate correlation functions that include these auxiliary fields explicitly, the 
extended BRST invariant stochastic action (5-5) is equivalent to the truncated one 
(5T1). Furthermore, the expectation values of local gauge invariant observables as 
functionals of V are independent of the stochastic gauge fixing functions and 0, 
and the F-P ghost loop effect gives no contribution to the renormalization in either 
(5-5) or (5T1), due to the retarded (or causal) property of their superpropagators. 
The decoupling of FP ghosts means that we may discard the FP ghost sector for the 
calculation of observables. 

Let us consider the generating functional of the stochastic Feynman diagrams 
for (5-5) by introducing the source terms 
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a 



(B-l) 



where Jq = Jd^zdr, = JdPzSr and J_ = JdPzdr. Here,^ denotes the appropriate 
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measure for the integrations with respect to the superspace coordinates and the 
stochastic-time. The dot ( " • ") denotes the sum with respect to the group indices. 
We also introduce the external source terms 

^rLj2l^..is^<^). (B-2) 

Here we note that 5brst</?" = (sf)", as defined in (5-4) and (5-7). The external 
source I^a has the same chirality as the original superfield (^". The transformational 
and dimensional properties of the fields and the external sources are listed in Table 
I. (In this table, the chirality + (— ) is defined by Dip = {Dip = 0).) 



Table I. Dimensional properties of fields and sources 
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The connected stochastic Feynman diagrams are generated by W{J, I) as 

ZiJ, I) = Jdfie''*''* ^ e^(-^'^) , i^tot = i^BRST + i^ex + Kf^^'^ , (B-3) 

where d/^ = 'DVVwVc'DcDc''Dc'V(pV(p'D BVB . The effective stochastic action, i.e., 
the generator of the 1-P-I vertices in the stochastic Feynman diagrams, is defined by 



r{cp,i) = w{j,i)-jY,J^'^^'" : 
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Hereafter, we omit "— "for the vacuum expectation values of the superfield (p°' = 
{V, m, c, c, c', c', (p, (j), B, B) so that these are not confused with their hermitian con- 
jugates. 

In terms of the effective stochastic action, the Ward-Takahashi identity for the 
extended BRST symmetry (5-4) is 

y:^^ = r*r = o. (B.5) 

^ Sif"^ 5I^c ^ ' 




In order to show the non-renormahzation theorem for the stochastic gauge fixing 
term, we derive local identities by changing the integration variable as c/ — c/ + (Jc/, 
with arbitrary 5c'. A similar identity is also derived for c' — c' + (5c'. In terms of F ^ 
these local identities are expressed as 

- * bl^ e 5Iv ) 5d 

— + -m^D^— - — = . B-6 
51^ 4r Sly j 5d ^ ' 

Similarly, iox B ^ B + SB and B ^ B + SB, we obtain 
Prom (B-7), the (S, S)-dependence of F is determined as 

where F does not depend on B, B, Ic' or 1^. The non-renormalization of the terms 
which depend on the Nakanishi-Lautrup fields in F is consistent with the power 
counting argument. Since F — F docs not include c', c', ly, or I^, F satisfies the 
local identities (B-6). By substituting (B-8) into (B-5), we obtain 



where (p" = {V,w,c,c,c',c',(f),(p). Let F^^iy and /r be the divergent part and the 
finite part of the effective stochastic action; i.e., /" = /r, + /div We assume that 
the loop expansion gives /div = Yle=i(^y^dtl determine F^^ from the Ward- 
Takahashi identity. In the following, we omit the superscript "(1) ". The local 
identity for /div indicates that /di^ takes the form 

p _ pi , pFP I p(cx) 
div — ^ div ' div ' div ' 
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where (sV") has the same dimensional properties as (s<y?") given in Table I. -T^;^ 
is a functional of (y, ro, (^). Here, in deriving (B-10), we have used the fact that 
the dependence of the divergent terms in Jjiv on the external soTirccs I.^a is at most 
linear. From the dimensional analysis, one might suspect that there exist terms 
quadratic with respect to, say, I^Sj'c'^ and I^TJ'c^, where Sj;' and TJ* are arbitrary 
functions of V. As we have already explained, because of the retarded (or causal) 
property of the superpropagator {Vw) with respect to the stochastic-time, a loop 
which consists of only {V w) as the internal lines does not contribute to the stochastic 
effective action. Therefore, the appearance of a term such as I^Sj'c^I^T^'^c'^ is 
excluded. Applying the same reasoning and the power counting argument, it can 
be shown that F^^^ is linear or quadratic with respect to the canonical conjugate 
momentum zu. In particular, there are no tu-indcpendent terms in I^j^ in (BTO). 
This is an essential point for the renormalizability of the SQM approach in terms of 
the BRST invariant framework with the stochastic action principle. 

The same conclusion can be derived from the so-called stochastic Ward-Takahashi 
identity. To derive this identity, we consider the generating functional Z (or W) in 
(B-3) under the condition Jy = = = /y = I^j, = = = 0. This condi- 
tion is required to perform the integration with respect to w without producing the 
non-trivial Jacobian factor. Integrating out vj by using the same trick as in (5T), 
we obtain 

Z\jy=j^=j^=iy=i^=i^=i^=o = independent of . (BTl) 

This means 

SG 

w = -— = 0, {or Jv = J^ = J^ = and Iv = I^ = U = 1^ = . (B-12) 

Therefore, in terms of the effective stochastic action, we obtain the stochastic Ward- 
Takahashi identity 

sr sr sr 

^ = ^ = ^ = 0, for tt7 = and /y = = = 7^ = . (B-13) 

This implies that there is no tu-indcpcndent term in i^jy- 

Let us decompose the Ward-Takahashi identity (B-9) into two parts. To do this, 
we assume the most general form of (sV") for F^f^ in (B-10). Prom the transfor- 
mational and the dimensional properties listed in Table I, the most general form of 
{s'ip"-) may be given by 

s'Wa = ^^U^lWcC^ - ^-W^lwJ' , 

s'(jf = 72kc" + Pg[(t) X c]" , 

= j2kc' + Pg[^ X c]" , (B-14) 
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where /3, (3 and 7 are divergent coefficients. The operators of dimension 0, namely 
^1 ~ -^"^b' ^ab ^^'^ ^ab ~ ^ oh ^ arbitrary functions of V. By using this 
definition, (B-9) yields two conditions. One is 

ss' + s's = Q. (B-15) 

The other is the condition to determine /div 1/^=0 = -^div + -^div ' 

/E{%^(»v") + %^(..")} = o. (B.16) 

where K\i^=q = i^BRST|B=B=o- Furthermore, (B-16) is decomposed into two parts, 

{V,zu,(t>,4>) 

E {^(»V") + ^gf(».")}=0, (B.17) 

{V,w,<t>,(l>,c,c) 

due to the (c', c')-dependencc. Here we comment on the generality of (B-14). Since 
and 0" are dimension-2 quantities, there is a possibility to exchange 4> (or cf)) 
and w in the definition (B-14). In the definition of s'cjf, the term g[zu x c]" is 
excluded because of the chirality of sV, i-e., D{s'(l)) = 0. Because of this chirality, 
this term cannot be multiplied by an arbitrary function of V. In the definition of 
s'tu", the possible terms might be U^lS'flf and W^^c?^"' where U^l and W^^ are 
dimensionless arbitrary functions of V. However, from the first condition in (B-17), 
since K includes terms linear with respect to ro, this implies that F'^^^ must include a 
Tu-indcpcndcnt term, and this is not the case, as we have explained above. Although 
this can be checked directly by adding these terms in the definition of s'ro" given in 
(B-14), to simplify the analysis, we discard this possibility in the definition (B-14). 

We first solve the condition (B-15). By substituting the definition (B-14) into 
(B-15), we obtain 

^ = /3 , 

X-^^d,L\^ + L\,d,X\^ - igf'\f5L\ + X« ) = , 

F|,5eL,«] + L^^d.Yl^ - igf^%f5L^^ + F^J = , 

Y'^d,L\ + Lid,X% - {X%dcL/ + L%dcY-a) = , (B-18) 

and 

-X^dddaL'^^ + daL%U,t^ + igPf^'daL'^, = , 
+iY%daL% - daL\W£) - {X%dadaL,' - daLj'U£) = . (B-19) 
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In (B-18) and (B-19), we have used the partial anti-symmetrization defined by 
[ab...cd] = ^{ab...cd — db...ca). Prom (B-18), the solution of (B-19) is given by 



daX% , W,^ = daY\ 



(B-20) 



By using the Maurer-Cartan-type equations for L""^ and L^" given in (A-6), we find 
the solution to (B-18) 

X\ = aL\ + {13- a)V'd,L\ , Y% = aL,'^ + {13- a)V''d,LC , (B-21) 

where a is a divergent coefficient. In terms of the renormalized BRST transforma- 
tion (5brst¥'" = A(s(/p"), (B-14) with (B-18) and (B-19) is expressed as 

s'V = a{sV) + {P- a)V''dc{sV) , 

s'w = a{sw) + {P- a){{sw) + Vddsw)} , 

s'c = P{sc) , s'c = (3{sc) , 

s'ct> = j{scj>) + {(3- j)g[cl) X c] , s'^ = 7(s^) + (/? - j)g[^ x c] . (B-22) 

Next, we solve the condition (B-16) and determine -Tdiv By substituting (B-22) into 
(B-17), it becomes 



+ 



+ 



(B-23) 



where s'V = V''dc{sV'') and = V''dc{szu''). The general solution of this 

equation consists of two parts. One is the non-trivial cohomology of the BRST 
transformation, i.e., the solution of '}2ae(Vzj 4> ip)"^ i^f"'){^-^div/^'^°') ~ ^- "^^^ other 
is a particular solution of (B-23). 

A particular solution of (B-23) is given by 



+ 



(B-24) 



This can be directly checked by using the commutation relations 



Jo SV- 



BRST 







, 



BRST 



(B-25) 
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where (^brst ~ ^aeiv-ua (p 4>)'§' {sip^)5 / dip"" . For the non-trivial cohomology of the 
BRST transformation, there exist three types of possible forms which are quadratic 
and linear with respect to w. We assume their most general forms, 

[e'''w''w'' , , fwa{n\V'' + Q\<t>'' + Q''l4>'') . (B-26) 

Jo Jo Jo 

Here i7"^ H", 0% and f2%, which are arbitrary functions of V, have dimensions 0, 2, 

and 0, respectively. It is not difficult to show that X""* cx C^ i?"^ oc S%, Q\ oc L\ 
and oc L^. In particular, in the third invariant, the relative numerical factors 
are fixed by the BRST invariancc. For S in the second invariant in (B-26), we obtain 

[d,L\E^ - L\dcS'')c' - (9,V^' - V^c^^Oc" = . (B-27) 

In order to clarify the meaning of this condition, we redefine S'^ as E°' = L'^jj^i and 
H« = Lj^H^. Then (B-27) reads 

c^ff = -[ c, Hi ] , c^ff = , (B-28) 

and 

c^ff = , c^ff = [ c, H2 ] , (B-29) 

where ) = ^^L^Ob, = j-yLj>db and S^2) = ^"(a)*"*- Now, (B-28) and 

(B-29) simply mean that Si and S2 are transformed as Si — > e'^^^'^Sie~'^^^'^ and 
S2 — e^s-^c^^g-igAc ^nfjep i]^Q BRST transformation of the vector superfield e^^^ — 
^ig\c^2gV ^-ig\c ^ Such Operators of dimension 3/2 are the so-called gluino fields, Wq 
and Wet- To construct the operators of dimension 2, one more local gauge covariant 
spinor derivative, = e~'^^^ D^e^^^ or = c^^^^qC""^^^, is necessary. There- 
fore, the second invariant is given by Hi oc {P", Wq} and S2 oc {^d, W^"}. Because 
of the reality condition, these two terms are not independent. The uniqueness of the 
operators Si and S2 of dimension 2 is confirmed as follows. The operators of di- 
mension 2, which do not depend on the stochastic-time derivative, are constructed 
from e^^^, e~^^^, and Da- They must include four spinor derivatives. For the 
transformation property, they also must consist of the same number of e^^^ and 
g-2gy particular, because of the chirality of the FP ghosts c and c, the spinor 
derivatives that exhibit the proper transformation property are T>a and D^^ for Si. 
Therefore, possible forms of Si can be constructed by successive operations of these 
derivatives on O, which is a quantity transformed as O ^ g*9-^c^g-tgAc gj^^gQ 
commutation relations of four spinor derivatives, using Va and Da, are expressed 
by Wa}, and V^Vm, where {P^, Da} = -2i(j™^P^, it is found that 

{P", Wa} is the unique possibility for Si. We can also conclude that {Pd, W"^} is 
the unique possibility for H2. 

In this way, we obtain the general solution of the reduced form of the Ward- 
Takahashi identity (B-23) 




32 



authors' name 



Here, K'^^\ and K^^^ are defined by 



k(i) = 



1 

WaWb , 

q2k 



K('^ = j^a{iV'' - ^(A/ - ^'L,-)} ■ (B-31) 

The divergent terms in the effective stochastic action are given by /div = -T^;^ + 
-^dhT^ -^dw'^' where the other two terms are defined in (B-10). The structure of 
the divergences 7^;^ is similar to the tree-level stochastic action, and thus we finally 
show that the divergences are absorbed by the renormalization of the wave functions 
and the coupling constants. 

Let us consider the scale transformation 

g e-^'Pg, k e^^^^^AC . (B-32) 
Under the condition 

X = ^{-a' + b' + c') , y = ^{a' - b' + c') , z = ^{-a' - b' + 3c') , (B-33) 

we have 

Then, using this relation, t'^^^ can be written as 

+(. + /3 - T)(/r^ + Z^"^) - + (. - .).§^}k (B.35) 

This expression can be expressed in terms of i^tot = -f^BRST + K^^^"^ . After some 
simple algebraic manipulations, we obtain 

f = f' J- f(FP) _i_ r(ex) 
-t div — div ^ div "T ^ div 
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d d d 1 

~^^dg ~ ^^'^Ok + (" ~ 7 + 2^ - ^)^^J -^tot • (B-36) 

This equation completely fixes the renormalization constants, which are introduced 
as 

^0 = l<p"0 = Zyjl^c^, go = Zgg, kq = Z^k, = Z^^, (B-37) 

where = {V,<f,c,c,c',c',(l),4>,B,B). By the definitions of the renormalization 
constants, (B-36) can be written 

^tot(^0' Vo>50,KO,Co) = Ktot{(p",I,p<',g,K,0 - ^div 

l-Zl/^ = x + P-a , l-Zl/^ = y , l-Zf = l-Zl/^ = z + P-^ , 

1 - zi/^ = 1 - 4/' = 7 - ^ - X = 1 - 4'^' = 1 - zy^ , 

1- Z^/^ = 1- Z^^ = P + x , 1- Z^^^ = a-j + z-x , 

ZipcZjipa = 1, for {V,w,c,c,c',(^,4>,^) . (B-38) 

In particular, we note that the relation Z^<xZi^a = 1 ensures that we can renor- 
malize the divergences of the effective stochastic action within this multiplicative 
renormalization procedure in the arbitrary order of the perturbation. 

The independent divergent coefficients are {x,y,z) ~ {a',b',c') and (a,/3, 7), 
while there are nine renormalization constants. This implies that there are three 
relations for the renormalization constants. In fact, we have 

Z-'Z'J^Z'J^ = 1, Z'J^ = Z^ZI/^ = Zy\ Zb = Zq = Ze = Z-e . (B-39) 

The first two relations are a consequence of the non-renormahzation of the stochastic 
gauge fixing term. The last one implies that the BRST transformation sd = iB 
{sc' = iB) is not renormalized. 

The analysis given in this appendix completes the proof of the renormalizability 
of the SQM approach in terms of the BRST invariant stochastic action principle in 
a formal sense. In the actual calculation, since the interaction is non-polynomial in 
the superfield formalism, we need to apply the background field method (BFM) . We 
note that, in BFM, the relations for the renormalization constants may be modified, 
for example, ZgZlf^ = 1 holds in BFM. 
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Although, we have not given a similar proof for the truncated version, (5-10) and 
(5-11), it could be possible to provide such a proof in terms of the Ward-Takahashi 
identity, which is also given by a form similar to (B-5). 

Iy^lL^ = r,r = o, (b-4o) 

where = (V, zu, c, c, c', c'). The external source Ii^a is coupled to sip°', which is de- 
fined by the truncated BRST transformation (5-10). In this case, the local identities 
for the connected stochastic Feynman diagrams, which are derived by the change 
of the integration variables (c' ^ c' + Sc' and c' ^ c' + 5c'), cannot be converted 
to those for the effective stochastic action. This indicates that the stochastic gauge 
fixing terms may be renormalized in the truncated version (5-11). If this is indeed 
the case, it would complicate the standard argument for the renormalizability. In 
particular, for the part of the stochastic gauge fixing and the F-P ghost sector, we 
would have to treat the most general form of the gauge fixing, because the truncated 
BRST symmetry is independent of the specified form of the stochastic gauge fixing 
functions. For example, in the renormalization procedure of the truncated version 
(5-11), it seems that there is no restriction on the modification of the gauge fixing 
functions, ^ = (p - {^D'^D'^h{V) and ^ = + {^D'^D'^h{V), where h(y) is an arbi- 
trary function of V. However, we emphasize that the BRST symmetry ensures that 
the expectation values of the local gauge invariant observables are independent of 
the stochastic gauge fixing functions and that the FP ghosts do not contribute to 
the renormalization. 
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